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I ' Abstract. We report first results of our systematic investigation of the 

magnetic field structure of rotating single magnetic white dwarfs and of 
white dwarfs in magnetic cataclysmic variables. The global magnetic field 
distributions on the isolated white dwarf HE 1045-0908 and the accreting 
^> ' white dwarfs in EE Eri and CP Tuc have been derived from phase-resolved 

. flux and polarization spectra obtained with EORSl at the ESO VLT using 

I the systematic method of Zeeman tomography. 
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1. Introduction 

Q . Magnetic fields have been detected in ~ 5 % of the 2300 known single white 

I dwarfs with field strengths ranging from 0.1-1000 MG and peaking around 

^ . 16 MG (McCook & Sion 1999, Wickramasinghe & Eerrario 2000). A similar frac- 

tion of magnetic white dwarfs has been found in accreting close binary systems 
(cataclysmic variables), possibly with a more restricted range of field strengths 
(7-230 MG, Beuermann 1998). The origin of the magnetic fields is not well 
understood. The magnetic fluxes, which should be conserved during stellar evo- 
\ lution, are similar to those of the magnetic main sequence stars such as the Ap 

and Bp stars which constitute ~ 5 % of the normal main sequence stars. The de- 
cay times of the lowest multipole components are predicted to be long compared 
to the evolutionary ages of the white dwarfs. The magnetic field topologies, 
at least of isolated white dwarfs, are, therefore, likely to be relics of previous 
evolutionary phases. In accreting systems, the field structure in the outer layers 
of the white dwarf, however, may have been significantly changed if the accre- 
tion rate is high enough that accretion occurs more rapidly than ohmic diffusion 
(Gumming 2002). 

Observational evidence suggests that the field topology of isolated and 
accreting white dwarfs can deviate drastically from that of a centered dipole 
(e.g. Wickramasinghe & Eerrario 2000, Schwope 1995). Information about the 
strength, orientation, and structure of the surface magnetic field of white dwarfs 
so far has been mainly derived from the analysis of photospheric Zeeman absorp- 
tion lines and thermally broadened cyclotron harmonics from the polar regions 
of accreting white dwarfs. The interpretation of Zeeman and cyclotron intensity 
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spectra alone is, however, often ambiguous and circular polarization spectra are 
required to investigate the magnetic field structure in detail. 



2. Zeeman tomography 

We have developed a systematic method, called Zeeman tomography, to derive 
the global magnetic field distribution on rotating white dwarfs (Euchner et al. 
2002). At present, it utilizes a database of 46800 sets of flux and circular po- 
larization Zeeman spectra of homogeneous magnetic white dwarf atmospheres 
calculated with a code developed by S. Jordan for a broad range of field strengths 
\B\ = 1-400 MG, effective atmospheric temperatures T = 8000-50000 K, and 9 
field directions relative to the line of sight. All spectra are calculated for a 
surface gravity log g = 8 and use a simple limb darkening law which is indepen- 
dent of wavelength. The magnetic field structure is approximated by a multipole 
expansion of the scalar magnetic potential, using spherical harmonics with co- 
efficients I = 1, ... ,5 and m = for the zonal and sectoral periodicity of the 
field. The database models are integrated for the B distribution of 900 elements 
of the visible white dwarf surface at a given rotational phase. A least-squares 
optimization code based on an evolutionary strategy is used to reconstruct the 
multipole parameters from a set of flux and polarization spectra obtained at 
different rotational phases of the white dwarf. 

The Zeeman tomography method has been thoroughly tested with the re- 
construction of various field geometries from synthetic spectra (Euchner et al. 
2002) . Here, we present first results of our application to real observational data. 



3. Observations 

We have obtained spin-phase resolved circular spcctropolarimetry of a sample of 
isolated and accreting magnetic white dwarfs using FORSl at the ESO VLT dur- 
ing four observing runs between May 1999 and December 2000. The instrument 
has been set up with a Wollaston prism and a rotatable quarter-wave retarder 
plate. Exposures have been alternately taken with retarder plate position an- 
gles (j) = —45° and (p = -|-45° in order to correct for instrumental polarization 
effects and linear-polarization cross-talk during the data reduction. Spectra of 
the target star and comparison stars in the field have been obtained simultane- 
ously by using the multi-object spectroscopy mode of FORS. This allowed us to 
derive individual correction functions for the atmospheric absorption losses in 
the target spectra and to check the data for remnant instrumental polarization. 
The target spectra cover the wavelength range 3800 - 7900 A at a spectral reso- 
lution A/ A A = 440. A signal-to- noise ratio S/N ~ 100 has been reached for the 
individual flux spectra at a spin-phase resolution A$ ~ 0.2. All raw data have 
been reduced using the context MOS of the ESO MIDAS package. 

All our spectropolarimetric observations of magnetic CVs have been ob- 
tained while the systems were in a low-state of accretion. Hence, the Zeeman- 
split photospheric hydrogen absorption lines were clearly visible as the emission 
of the white dwarf was not diluted by emission from the accretion stream. 
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Figure 1. Flux (top) and circular polarization (bottom) spectra of 
HE 1045-0908 at (the arbitrary) rotational phases $ = 0.0, 0.25, and 
0.5. The synthetic spectra for the best-fit model (thick line) consisting 
of a dipole, quadrupole, and octopole field component is shown super- 
imposed on the observed spectra (grey curve). For clarity, the upper 
two flux spectra have been offset by 2 and 4 flux units, respectively, 
and the polarization spectra by 0.1 and 0.2 units, respectively. 



4. Results 

4.1. HE 1045-0908 

HE 1045-0908 is an isolated hydrogen-line white dwarf with a polar field strength 

B ^ 20 MG and a rotational period Prot ^ 2 — 4hr (Reimers et al. 1994, Schmidt 
et al. 2001). Our observations cover ~ 2hr during which drastic variability in 
the flux and polarization spectra has been seen (Fig. 1). The field topology 
can be described by a combination of aligned dipole, quadrupole, and octopole 
terms with a common offset perpendicular to the magnetic axis (Fig. 2). An 
almost equally well fit has been obtained for a multipole expansion including 
the zonal components up to / = 5. It is reassuring that both fits yield similar 
ranges for the observed field strengths and viewing angles. 
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Figure 2. Distribution of the total magnetic field strength B and 
its longitudinal component Bi on the white dwarf in HE 1045-0908 at 
(arbitrary) rotational phases $ =0.0, 0.25, and 0.5 (left) and frequency 
distribution of B and the viewing direction cosine cos$ (right). 



4.2. EFEri 

Spectropolarimetric observations of EF Eri have been obtained during two epochs, 
in December 1999 and in November 2000 and cover in total 2 spin cycles of the 
white dwarf. Rotational phases have been calculated using the ephemeris given 
by Piirola et al. (1987), with $ = 0.0 corresponding to the narrow minimum in 
the infrared light curve of EF Eri, i.e. to the phase when the magnetic pole is 
closest to the observer. As our data show little variation in the Zeeman features 
modeling has been done for the phase averaged flux and polarization spectra 
(Fig. 3 and 4). The data are sufficiently well described by a multipole expan- 
sion including the zonal components up to Z = 5. We note, however, that there 
are still systematic residuals in our best fit solution which indicate that the field 
topology must be even more complex. 

4.3. CPTuc 

CPTuc (= AXJ2315-592) has been discovered by Thomas & Reinsch (1996) 
who derived an upper limit of the magnetic field strength at the accreting pole 
B < 17 MG from the properties of the optically thin cyclotron emission. Our 
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Figure 3. Average flux (top) and circular polarization (bottom) spec- 
trum of EF Eri. The synthetic spectrum for the best-flt model (thick 
line) consisting of a five component multipole expansion is shown su- 
perimposed on the observed spectrum ( grey curve ). 

spectropolarimetric observations cover ^^0.5 spin periods of CPTuc (Fig. 5 and 
6). Rotational phases have been calculated using the ephcmcris given by Ramsay 
et al. (1999) where = 0.0 is defined by the minimum of the X-ray light curve. 
Our best-fit solution for the field topology comprises a dipole, quadrupole plus 
octopole combination. Again, residual features indicate an even more complex 
field. The high polar field strength found here is not contradicting the much 
lower upper limit derived by Thomas & Reinsch (1996) for the field strength in 
the accretion region as the latter is probably offset from the magnetic poles. 



5. Conclusions 

We have shown that Zeeman tomography is a suitable systematic method to 
derive the global magnetic field distribution on rotating white dwarfs. With 
this method we have obtained for the first time detailed information about the 
range of field strengths and the field topology of a sample of isolated and ac- 
creting white dwarfs. Our results clearly demonstrate that a single value as 
obtained with hitherto available methods is not sufficient to quantify the field of 
magnetic white dwarfs. Overall our model fits are in excellent agreement with 



6 



K. Reinsch et al. 



filling factor: n^^^^^^™ 0.067 




Figure 4. Distribution of the total magnetic field strength B and its 
longitudinal component Bi on the white dwarf in EF Eri at spin phase 
<I> = 0.0 (left) and frequency distribution of B and the viewing direction 
cosine cos^> (right). 



observations. Remaining differences indicate that the field topology is even more 
complex than described by an up to 5 component multipole expansion. A more 
detailed discussion of our results will be presented elsewhere (Euchner et al., in 
prep.). 
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Figure 5. Flux (top) and circular polarization (bottom) spectra of 
CPTuc at spin phases ^ = 0.99, 0.21, and 0.43. The synthetic spectra 
for the best-fit model (thick line) consisting of a dipole, quadrupole, 
and octopole field component is shown superimposed on the observed 
spectra (grey curve). The upper two flux and polarization spectra have 
been offset by 0.2 and 0.4 units, respectively. 
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Figure 6. Distribution of the total magnetic field strength B and its 
longitudinal component Bi on the white dwarf in CP Tuc at rotational 
phases $ =0.99, 0.21, and 0.43 (left) and frequency distribution of B 
and the viewing direction cosine cos^ (right). 



